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ABSTRACT 


The  transient  swan-square  response  of  a nonlinear  s ingle- degree- 
of  free don  Mechanical  system  to  nonstationary  random  excitation  char- 
acterized by  the  product  of  an  envelope  function  and  a stationary 
Gaussian  random  process  is  determined  by  the  equivalent  linearization 
technique.  A unit  step  envelope  function  is  considered  in  conjunction 
with  both  correlated  and  white  noise  with  zero  mean. 

It  has  been  shown  that  for  white  noise  modulated  by  a unit  step 
function,  the  transient  mean-square  response  never  exceeds  the  stationary 
response.  However,  the  mean-square  response  to  correlated  noise 
aodulated  by  a unit  step  function  may  exceed  its  stationary  value. 

Thn  analysis  is  extended  to  the  multi -degree-of- freedom  nonlinear 
system  for  the  case  of  mutually  uncorrelated  noise. 


CHAPTER  I 


INTRODUCTION 

The  transient  naan- square  response  of  a linear  slngle-degrte-of 
fre^dcm  Mechanical  sys*«M  to  certain  types  of  nonstationary  random 
excitation  has  been  studied  by  several  authors  [1,  2,  3,  4].  The  non- 
stationary Input  was  taken  In  the  fo>??  of  a product  of  a w*ll -defined 
envelope  function,  A(t)  and  a stationary  Gaussian  noise  with  zero  Mean, 
n(t). 

T.  X.  Caoghey  and  H.  <3.  StuMpf  [1]  have  exaadned  the  case  In  «nlch 
tbe  envelope  function  A{t)  was  a unit  step  function  and  n{t)  was  assumed 
to  be  either  white  noise  or  f;roadi-band  noise  whose  power  spectral  desn- 
sity  has  no  sharp  peaks.  Results  of  their  analysis  were  applied  to  the 
<tete miration  of  the  structural  response  to  earthquake  ground  Motion. 

V.  V.  B-:  ii1n  [2]  hiss  determined  the  Mean-square  response  of  a linear 
structure  represented  by  a second  order  differential  equation  when  shat 
structure  is  subject  to  earthquake  axel  cation.  In  his  analysU,  he 
considered  the  ground  acceleration  to  be  characterized  by  the  product 
of  an  exponentially  decaying  ha  neon 1c  correlation  function  and  an 
envelope  function,  A(t)  * Vc*. 

In  a recent  paper  [3],  R.  L.  Bamoski  and  «J.  R.  Maurer  have  forwu- 
lated  the  tine  varying  ns^-square  response  of  a linear  single- degrte- 
of-freedo«  system  In  terns  of  the  system  frequency  response  function 
and  the  generalized  spectral  density  function  of  the  input  excitation. 
They  cons  1C  red  the  envelope  function  to  bo  either  the  unit  step  function 


i m Jill  L 1 1 .y.-^Tfr  > w-*ri  > 


m 


or  a rectangular  stop  function.  L.  L.  Bucciarelli  and  C.  Kuo  [4]  bum 
recently  obtained  an  approx  inete  across  Ion  far  tbe  ■can  square  rat  paw 
to  oxcilatiw  characterized  by  a general  envelope  function  subject  only 
to  tbo  restriction  ithat  tbe  envelope  function  Is  slowly  varying.  Their 
work  also  pave  an  astiaafad  maxim*  value  of  tbe  tan  square  response. 

In  all  tbe  above  studies,  tbe  systavs  treated  were  linear. 

To  date,  tbe  p rob  Ian  of  response  of  a nonlinear  systen  to  nen- 
stationary  random  excitation  has  beat;  tntiened  in  only  one  piece. 

There,  It.  H.  Toland,  C.  Y.  Yang,  and  €.  Hsu  [5]  anployed  a randan  walk 
nodal  to  determine  systen  response  to  stationary  touts  1 an  white  noise. 

t 

Tbo  extension  to  tbe  nonstatlomry  case  was  discussed,  but  not  carried 
through  to  oonpletion.  There  ere  nany  sys tarns  whose  notions  are 
characterized  by  nmrlzmr  differential  equations,  particularly  ufcen 
tbe  notions  ore  large.  It  Is  tbe  purpose  of  this  study  to  presont  an 
approxlnate  solution  to  tbe  transient  naan- square  response  of  a staple 
noel  iwer  systen  to  e nonstationary  randon  excitation.  Only  systens 
with  gtowtric  nonlinearities  (rather  than  naterials  nonlinearities) 
Involved  are  considered  end  tbe  nonlinear  differential  equation  is 
linearized  by  an  equivalent  linearization  technique.  All  conputetion 
for  obtaining  the  wan- square  response  of  tbe  system  is  straightforward. 
However,  it  Is  not  generally  possible  to  obtain  a closed  form  solution 
of  tbe  neon-square  response  and  tbo  aquation  has  to  bo  solved  numerically. 
First  a $1ngle-degr*i^f- freedom  system  is  treated  In  Chapter  II  and  then 
a multi  d*prw-of-f>~*uga  system  Is  discussed  for  various  special  cases. 

Results  of  he  present  Investigation  could,  as  an  approximation, 
be  applied  to  response  determination  of  thin  elastic  plates  and  shells 


[3] 


(regarded  as  single  degree  of  freedo*  system  for  any  given  node) 
when  these  system  are  subject  to  pressure  fields  that  excite  large 
amplitude  oscillations.  Jet  engine  sound  pressures  would  be  one 
example  of  such  a pressure  field. 


[4] 
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CHAPTER  II 

A SINGU-DEGREE-QF-FREEDOM  NONLINEAR  SYSTEM 


2.1  Statement  of  the  Problem 

Consider  a lightly  damped  si ngle-degree-of- freedom  mechanical 
system  subjected  to  a random  excitation  and  governed  by  the  equation 
y(t)  + 2^{t)  + w2n(y(t)  + g(y) ) » f(t)  (2-1) 


where 

C * fraction  of  critical  damping 

u>n  * natural  frequency  of  the  corresponding  linear  system 


g(y) 


H J2k  + 1 

T. 


uk  > 0 


k«l 

The  nonstationary  random  excitation  f(t)  is  expressed  by 


(2-2) 


f(t)  - A(t)n(t) 


(2-3) 


where  A( t)  is  a well-defined  envelope  function  and  n(t)  is  a Gaussian 
stationary  random  process  with  zero  mean  and  autocorrelation  function 

Rn'T>' 

O 

We  are  to  determine  the  mean-square  response  E[y  (t)]  to  an 
input  f(t)  when  the  envelope  functions  are  a unit  step  function  and 
exponential  function,  respectively,  i.e.. 


A( t)  * u(t)  (2-4) 

M 

A(t)  » Z A.exp(-c,t)u(t) 

1-1  1 1 


(2-5) 


[5] 


and  n(t)  has  tl*  following  autocorrelation  functions 


Rn(T)  * 2itk06(i) 

(2-6) 

Rr(t)  = K0exp(-aj-rj)cosBT 

(2-7) 

where  u(t)  Is  a unit  step  function  and  5(t)  is  the  Dirac  delta  function. 
Note  that  If  we  let  all  A^  be  zero  except  for  A^  « 1 and  c^  * 0,  Eq.  (2-5) 
reduces  to  Eq.  (2-4). 

2.2  Response  Formulation 

Although  various  methods  can  be  applied  to  determine  the  response 
of  nonlinear  systems,  the  equivalent  linearization  technique  will  be 
used  here.  This  technique  was  developed  by  Krylov  and  Bogoliuvov 
for  the  treatment  of  nonlinear  systems  under  deterministic  excitations, 
and  then  R.  C.  Booton  [6]  and  T.  K.  Caughey  [7]  applied  this  technique 
to  problems  of  random  vibrations. 

He  assume  that  an  approximate  solution  to  Eq.  (2-1)  can  be  obtained 
from  the  linearized  equation 

y + 2^  + w^y  * f(t)  (2-8) 

2 

where  Be  Is  the  equivalent  linear  damping  coefficient  and  u»c  Is  the 
equivalent  linear  stiffness.  The  error  "e"  due  to  linearization  is 
given  by  the  difference  between  Eqs,  (2-1)  and  (2-8),  l.e. , 

e * 2(cwn  - 8e)y  + (wjj  - u^Jy  + g(y)wjj  (2-S) 

2 

The  variables  Be  and  are  selected  so  as  to  make  the  mean-square 

2 2 
error  E[e  ] a minimum.  The  minimization  of  £[e  ] require  that: 


[6] 


5^1  * 0 (2-10) 

58e 

3E[el  , 0 

Substitution  of  Eq.  (2-9)  into  Eqs.  (2-10)  and  (2-11)  gives 


2 

• -»(wn-ee)E[y2)  - 4W2  - <.2)E[yy] 

- «Cyg(y)3u>*  - 0 (2-12) 


2 

* "4(?w  “8  )E[yy]  - 2{<l  - o>2)E[y2] 

aUp)  n 6 

- 2E[yg(y)]u>2  «=  0 (2-13) 

2 

Solving  for  8e  and  oie  from  Eqs.  (2-12)  and  (2-13),  we  have 


a = 25.  ♦ u2  EEyz3ELya(y)3  - SLaJilaMJ 

e " n E[/]E[}2]  - (Ebiir 

u2  . + „2  E[y2]Eto(y)]  - E[yy]E[yg(y)] 

e " n E[yr]E[>2]  - (E[yj])* 


(2-14) 


(2-15) 


From  Eqs.  (2-12)  and  (2-13), 


(2-18) 


= 4ECyy] 

3Be3(wg) 
and 

ifllfiL  sW| . = l6(E[y2]E[y2]  - (E[yy])2} 

36g  8(0^  36e3(u)g) 

= 16det(K)  (2-19) 

where  det(K)  Is  the  determinant  of  the  matrix  of  covariances.  Since 
the  upper  bound  for  the  nonstationary  cross-covariance  E[yy]  Is  given 
by  the  Inequality  [8], 

(E[yy])2  < Ety2]E[y2]  (2-20) 


then. 


&[e2] 
SBe  W‘) 


2 


From  Eqs.  (2-16),  (2-17),  and  (2-21),  It  can  be  seen  that  the  conditions 
(2-14)  and  (2-15)  truly  give  a minimum  E[e2]. 

In  order  to  express  the  right  hand  side  of  Eqs.  (2-14)  and  (2-15) 

In  terms  of  E[y2],  E[y2]  and  E[yy]»  It  Is  necessary  to  know  the  prob- 
ability density  function  p (y,y).  In  general,  hotever,  p(y,y)  Is  not 
known.  If  the  Input  Is  Gaussian  and  the  nonl Inearl  ties  of  the  system 
are  small,  then  the  response  of  the  linearized  equation  (2-8)  Is  also 
Gaussian.  Therefore,  the  assumption  Is  made  that  the  probability  density 
function  p(y,y)  Is  Gaussian  with  covariances  to  be  determined. 


08] 


Before  constructing  the  probability  density  function,  the 

enseufcle  average  of  y and  y is  calculated  by  use  of  Buhners  Integral. 

Assuming  that  the  system  is  at  rest  initially,  we  have  the  solution  of 

Eq.  {2-8)  to  be 
t 

y{ t)  * j h(t-T)f(t)dr  (2-22) 

0 

Where  h(t)  is  the  impulse  response  of  the  system  defined  by 

,-v 

h(t)  - (Slikidt)u(t)  (2-23) 

d 

and 

wd  * we  “ $e  (2-24) 


The  ensemble  average  of  y is  obtained  by  taking  the  ensertrie  average 
of  Eq.  (2-22). 

E[y]  « ^ h(t-T)E[f(t)]dr  (2-25) 

Since  we  assumed  that  E[f(t)]  * 0,  then 

E[y]  « 0 (2-26) 

Similarly,  the  ensetfcle  average  of  y is  obtained  by 


E[y] 


^ h{t-T)E[f(T)]dr  - 0 

' ) 


(2-27) 


Because  of  Eqs.  (2-26)  and  (2-27),  the  assumed  Gaussian  probability 
density  function  p(y,y)  takes  the  form: 


W— •'*  -'.^s  ,***►->-  -v. , 


I 

t 


^*«wt"»  -^■wrH**<^«*^t>*s-  - 


winr<«  v *s« 


[9] 


p(y,y)  * - r exp(-ay2  + 2byy  - cy2) 

2ir(det(K))* 

where 

a * E[y2]/(2det(K)) 

b * E[yy]/(2det(K))  . 

c * E[y2]/(2det(K)) 

det(K)  - E[y2]E[y2]  - (E [yy])2 


(2-28) 


(2-29) 


Expression  for  E[yg(y)3  and  E[yg(y) ] are  obtained  In  the  Appendix  and 
r.re  shown  to  be  [9], 


E[yg(y)3 


ys(y)p(y.}/dydy 

< 

mC9 


• ? ^ (2-30) 

k-1  ^ rkl 

I 


It  1$  interesting  to  observe  that  the  above  equivalent  linear  diaping 
2 

26«  and  stiffness  w.  are  Identical  to  those  found  for  a stationary 
process  In  which  case  E[yy]  Is  equal  to  zero  and  Eqs.  (2-14)  and  (2-15) 
are  simplified  as 


2Wn  ♦ 


Eialill 

eS*T 


UaSiLl 

Ely8] 


(2-3*) 


If  Instead  the  nonlinearity  Is  Involved  only  In  the  velocity  term  such 
as  <j(y)»  we  can  easily  show  that  the  equivalent  linear  damping  and  stiff- 
ness for  a nonstationary  process  are  the  sane  as  those  for  a stationary 


process. 

2 

The  nean-square  response  E[y  ] at  any  Instant  of  tine,  t.  Is 

2 

obtained  fron  the  computation  of  the  expected  value  of  (y(t)}  over 
the  enseafcle  of  the  response.  Fron  Eq.  (2-22), 


Ely2]  ■ 


ftft 

h(t-T)h(t-T)E[f(T)f(t)]drdr 

0 Jo 


(2-35) 


Since  A(t)  is  well-defined  function  and  n(t)  is  stationary. 


Etf(T)f(x)]  - A(T)A(t)E[n(x)n(T)] 

- A(t)A(t)Kh(t-t)  (2-36) 

where  Rn(T-T)  Is  the  autocorrelation  function  of  n(t).  Substitution 
of  Eqs.  (2-23),  (2-32)  and  (2-36)  Into  Eq.  (2-35)  leads  to 


[11] 


E[y2]  ■ 


t ft 


h{t-T)h(t-t)E[f{T}fCT)3dTdT 


0 Q 


l I 


X sin*»d(t-T)A(T)A{T)Rn(T-T)drdT 


Fran  Eqs.  (2-24),  (2-32)0  and  (2-33),  wj  becomes 


(2-37) 


. 1 ♦ I il.  13UJ1S.  <E[y?])k 
k-1  ^ 2*k! 


(2-38) 


With  RJt)  defined  by  Eqs,  (2-5)  or  (2-7),  a»0  A(t)  defined  by  Eqs.  (2-4) 

n 

or  (?-5}0  Eqs.  (2-37)  end  (2-38)  becos«>  sioultaneous  nonlinear  algebraic 
equations  for  £[y  ]. 

Alternatively,  Eq.  (2-37)  can  be  expressed  In  tarns  of  th*  power 
spectral  density  #(w)  of  n(t).  Since  the  autocorrelation  function 
VT  -?)  is  given  by 


»•* 

I 

^(;«j)coa«(T-T)dw 


(2-39) 


</  •' 


Then,  substituting  Eq.  (2-39)  Into  Eq.  (2-37)  and  hanging  t ho  order 


of  Integration,  we  haw 


Ety2] 


* ft  ft 


rm  t*r 

jUO 

0 * J.  .. 


«xp{-cu>n(  2t-T-t  }s1nw^(  t-x) 


x s1nh»j(t-T)A(t)A(T)cosw(T-t)drdtdu 


(2-40) 


If  the  Input  Is  asswned  to  be  white  noise,  then  substitution  of 
Eq.  (2-S)  lute  Eq.  (2-37)  leads  to 


E[y2l  * ^jr  [<«q>  -2c»||(t-T))*2(t)i1ii?M(|(t-T)<h  (2-41) 


2.3.1  Unit  Step  Envelope  Function 

Let  us  first  consider  the  case  In  which  the  envelope  function  Is 
defined  by  Eq.  (2-4).  By  perforating  the  Integration  of  iEq.  (2-41),  we 
obtain 


.2-t  . "*o 


20  3 ■ -f  txp{-2cwn(t-r)}sinfw|j(t-T)dT 


-2cwut  2c  wl  « 

jr_  [1-e  (1  ♦ — n — sin  Mjjt 

«*Si‘A£  ♦ “d>  -d  d 


♦ Jlnj.jt)] 


(2-42) 


[133 


Employing  Eq.  (2-38),  Eq.  (2-42)  becomes  a nonlinear  algebraic 
equation  for  E[y2]  since  is  a function  of  E[y2]  In  Eq.  (2-38).  This 
type  of  equation  generally  has  more  than  one  solution.  However,  from 
physical  considerations  the  cksfred  solution  will  be  that  one  close  to 
the  solution  of  the  corresponding  linear  system  because  inly  a weakly 
nonlinear  system  Is  betag  considered.  Since  the  general  procedure  for 
solving  the  nonlinear  algebraic  equation  Is  not  available,  we  shall 
use  Newton's  method  of  tangents  to  obtain  an  approximate  solution  at 
Instantaneous  values  of  tine,  t,  and  then  Iterate.  The  solution  by 
Newton's  method  sometimes  does  not  converge  If  a poor  Initial  value  Is 
cx>sen.  However,  since  the  mean-square  response  of  the  linear  system 
E[y03  Is  assumed  to  be  close  to  that  of  the  nonlinear  system,  E[yjj3  1? 
suitable  for  use  as  the  Initial  trial  solution  for  an  Iteration  scheme. 
Throughout  the  present  study,  this  Iteration  scheme  together  with 
Newton's  method  Is  used  for  obtaining  an  approximate  solution. 

As  a numerical  example,  let  us  consider  the  simplest  case,  l.e., 
g(y)  « uy3  (2-43) 

For  various  values  of  y end  damping  coefficient  c,  E[y  ] is 
computed  and  the  normalized  plots  are  shown  In  Figure  1 through  3. 

The  normalization  factor  Is  determined  by  the  stationary  mean-square 
response  of  linear  system 

EEyfA  ■ (2-43) 

The  parameter  y Is  chosen  in  such  a manner  that  given  y the  stationary 
mean- soo* re  response  reaches  40  percent 0 60  percent  and  80  percent  of 


[14] 


2 

E[yPs«  If  the  damping  Is  small,  Eq.  (2-42)  can  be  approximated  by 


1 - e 


-2suit 


E[y2]  * E[y2] ty 

0 s i + 3pE[/] 

from  which  the  following  approximate  solution  is  obtained. 

E[y2]*  ^ {[1  + 12uE[y2]s(l  - a’2^)]^!} 


(2-4S) 


(2-46) 


In  what  follows,  let  us  show  that  the  transient  mean-square 
response  for  both  linear  and  nonlinear  systems  does  not  exceed  the 
stationary  mean-square  response  to  white  noise.  That  is. 


e[y23«  i E[y2] 


os 

ECv2],  > E[y2] 
From  Eq.  (2-42),  we  have 

wK 


(2-47) 


E[y2]  * 


4^nwe 


7 


l-e^n^M^)2  + (-CJI)2+  (^H)^ 

w A Klj  'll)/ 


x sf n(2ti»rft  - 6)] 


(2-48) 


where 


0 - tan'1  (2^) 
wd 


Since 


1 - (ft*  > eft2  ♦ ft¥ 


u 


co„ 


CO. 


-1  < sin(2wdt-9)  < 1 


1 V2cwntCl+(^)2  ♦ y^)Z  + (^4  *1«(2»dt-0)]  < 1 (2-49) 


Tt»  equality  holds  for  t •*  «. 


Therefore, 

E[v2)  £ eK0/4^0+3^Ly2]) 
Solving  for  CCy2l»  It  Is  concluded  thet 


(2*50) 


wK 

E[y2]  -1  - E[y2]s  m 

l *CU|'  ) 


For  linear  systems,  substitution  of  p * 0 Into  Eq.  (2-50)  leads  to  the 
first  equation  of  (2-47). 


2.3.2  Exponential  Envelope  Function 

If  a white  noise  Is  nodulated  by  the  exponential  envelope  function 
described  by  Eq.  (2-5),  then  Eq.  (2-41)  becomes 


E[y2]  - 


nr 

wd  'o 


exp[-2&»n(t-x)]sin  «d(t-t) 


a M 

X I J AjAj8Xp(-(Cih!j)T]dt 

i«rf  1 J 1 J 


-2tw_t 

N M KJLA.e  n 

E E &-* y 

1*1  j*1  1 j"*t**d^ 


[e2Xijt-(l+— Ji*in2«dt 


(2-51) 


m 


♦ Sln2  .t)3 
wd  d 


where 

x1j  * -*(c1  * CJ> 


(2-51) 

cont'd 

(2-52) 


Consider  the  special  case  In  which  the  envelope  function  and 
nonlinear  tern  g(y)  are  the  following: 

A(t)  * Ae~ct 


g(y)  * hv3 


(2-53) 


Then,  Eq.  (2-51)  Is  simplified  to  the  fonw 


try2] 


S 


’"V 


e-2« 

-SB W 

« -2c«+x* 


,2(?-5)t.tH  2j^i  $1n2(l2.,2)i 
x -6 


where 


2 

x 

T 


6 


* x $1n2(x2-c2)*T] 
(xZ-«2)* 


I ♦ 3wE[y2] 

V 

c/“n 


(2-54) 


(2-55) 


Equation  (2-54)  was  solved  by  a numerical  Iteration  method  described 
In  Chapter  2.3.1,  and  normalized  plots  are  shown  In  Figure  4 and  5. 

The  normalization  factor  Is  chosen  In  this  caso  to  be 

4 m (2-M) 

which  is  the  stationary  naan-square  response  for  the  case  of  c « ft 
and  u » 0.  Two  cases  6 * 2t  and  6 • ft  ant  Illustrated  In  those  figures 
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and  the  results  show  that  for  large  6,  l.e.,  rapidly  decreasing 
amplitude  of  the  Input,  the  effect  of  the  parameter  u Is  Insignificant. 
If  c.  4 ♦ !*  Eg.  (2-54)  can  be  expressed  approximately  by 


e|>23  - a2 

Solving  for  E[y2] 

ECy2]  • l 

Is  obtained. 


(t-«)0+3«E[yS3) 

{[H^O  (e-2«^-2cT);ji  _1} 


(2-57) 


2.4  Response  to  Shaped  Correlated  Noise 

If  the  Input  noise  Is  assumed  to  be  correlated  as  the  damped 
harmonic  form  of  Eg.  (2-7),  then  Eg.  (2-37)  becomes 


Efy2] 


*0,-2<V 

— 1 — 
“d 


rL  t 

t a • 

exp[cwn(T+r)“CtjT-T|  ]A(t)A(t) 


0 0 


• I I 

s 1 nw^  ( t-T ) s 1 nuj ( t- t ) cosg ( t-t ) dTdr 


(2-58) 


2.4.1  Unit  Step  Envelope  Function 

Substituting  Eg.  (2-4)  into  Eg.  (2-58),  we  find  that  Eq,  (2-58) 
takes  the  form  , 


ft 

> 

j 

Efy23  • -*-1 — 

exp[(cu»n-o)T]sinw(j(t-t)  - 

4 

J 

0 

. 4 

rr 


exp[(5wn+a)x] 


X SilSJ^t-TjcOfiBCT-TjdT 


ft 


dr  + 


■ « 

exo[  ( t<un+ah]$  in^(  t-t) 


j 


D 


[13] 


x < 


* 

exp[{cwn-a)T]s1nwe|(t-T)cos$(T-T)dT 


dr 


(2-59) 


The  double  Integral  in  Eq.  (2-59)  may  be  evaluated  after  some  tedious 
algebra  to  give. 


E[y2]  * -f  e‘2jV  l (UtVUt) 
-d  w 


(2-60) 


where 


Q1  « eplt(R1-R2)p1/2R1R2 


V 

q3  * 
Qs  * 

V 


eplt(ft1R24ft2R1)/2R1R2 
pt(r1+R2)  p2(R3+R4) 


TK^T 


Q1R4^2R3 

~WT 


q1r2+o2s1 

~*w 


P2 

q7  * * 

^2 

*>»  * ’ 3^ 


(2-61 ) 


1193 


T-j  = [p2ep2t(Ri|-R3)+R^(p2COsnit-fl^s1nn^t5  + 
R3{fl2s1nft2t-p2cosn2t)]/2R2R4 

I2  * [eP2t(n-jR^-K]2R3)-R4(p2sirtfi-jt+Q.jCOsn-jt)  - 

R3( p2si nn2t+o,cosn2t) ]/2R3R^ 


T3  * (1-ep3t)/2p3  + P3ep3t/2R5  - 
P3cos2uj£,t-2u(;|Sln2(Dcjt)/2Rg 
P3t 

-p3s1n2w£.t-2uj(jCOs2a)(jt)/2Rg 


Tg  * -epl*(p-jSlnOj  +p.|COsn.|t)/2R.j+p1/2R.j  - 

(n2s1nfi,  t-p1cosn1t)/2R2-(n2s1n2a)cft-p1cos2wc|t)/2R2 


Tg  * ep1t(ft1cosft1t-p1s1nfl1t)/2Rrn1/2R1  + 

eh^PiSinOjt+flgCos  -jt)/2R2  - 
(p-jSi  n2u(jt4fl2cos2uc|t)/2R2 

T?  * ep1t[(n1R2+Q2R1)sinn2ttp1(R1-R2)cosn2t]/2R1R2  ” 
[n1  R2S 1 n2todt-p  ^ ( R2cos  2wdt-  ) ]/2R^  R2 


(2-61) 

cont'd 
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i 

I 

i 

j 


« 

i 

t 


Tr  * ep1t(pxs1nn2t+n1cosfi2t)/2R1-n2/2R2  - 

(p^cos2u)dt+n-jCOs2w(jt)/2R^  - 

e^Cp-jSlnflgt  - 2cosw(jt)2R2 

pl  = *Va 
p2  “ 

p3  * p,  ♦ p2 

nl  * ud"^ 
fi2  * wd+8 

2 2 

R,  * P]  + 

2 2 

P2  ” P1  + ^2 

R3  * p2  + P1 

R4  * p2  + P2 

R5  * (Vp2*2-Mwd 

From  Eqs.  (2-60)  and  (2-61)  It  Is  seen  that  the  mean-square 

response  depends  upon  an  interrelationship  which  involves  the  system 

damping  c,  the  corresponding  linear  system  natural  circular  frequency 

wn,  the  decaying  constant  a and  the  frequency  8 of  the  correlation  functict 

2 

Determining  the  solution  E[y  ] requires  much  algebra  even  for  the  simplest 
case  of  g(y)  e My0. 


(2-61) 

cont'd 
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When  the  Input  Is  white  noise*  only  the  value  of  damping  of  the 
system  affects  how  quickly  statloaarlty  Is  attained  as  seen  In 
Section  2.4.  However,  for  a correlated  noise  Input,  the  time  required 
for  the  resoonse  to  reach  a stationary  value  Is  Influenced  not  only 
by  the  system  damping  coefficient  ; but  also  by  the  decay  constant  a 
of  the  Input  noise.  Normalized  plots  for  g(y)  yfcy  are  shown  In 
Flqures  6 through  8 for  various  value  of  a with  fixed  8.  The  normal 1- 
zatlon  factor  was  chosert  as  the  stationary  mean-square  response  E[y03s 
of  the  linear  system  which  can  be  obtained  by  letting  t » and  y * 0 
In  Eq.  (2-60).  Instead  of  finding  the  stationary  mean-square  value 
from  Eq.  (2-60),  we  obtain  It  from 


(*♦(0))  - 

^ II 


(wjj-w  )+l2yun<i>y 


(2-62) 


where  <fr(w)  Is  the  one-sided  power  spectral  density  of  n(t): 


(2-63) 


By  contour  Integration,  E[y*l  becomes 

s 


Er»2l  . S°°'  V1'  > W + . , 

0S  C(l-tZ)i(A|+»|) 


(2-64) 
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where 

Aj  « o,2+B,2+1-2c2+2B'(1-C2)* 

B,  • 2c[(l-c2)i+e‘] 

a2  « a,2+e,2*i-2c2~2B'(i-c2)* 

B2  - 2cCO-;2)i-B'] 

C « H(ft,2-e'2)2-4a,2B,2+2($,2-«,2)(2c2-l) 

D « 4a'Bs(1-2c2+a‘2-B'2)  {2-65) 

a*  a/u>n 

B'  B/a>n 

This  was  done  for  checking  purposes.  Since  the  normalized  mean-square 
value  must  asymptotically  approach  unity  for  large  t,  we  can  check  the 
results  of  Eq.  (2-00)  If  we  let  p*0.  Results  show  that  as  a decreases, 
l.e.,  the  power  spectral  density  has  a sharp  peak  at  some  frequency, 
then  the  transient  response  tends  to  exceed  the  stations*?  value. 
Another  Interesting  result  Is  that  the  nonlinear  response  becomes 
greater  than  the  corresponding  linear  response  under  certain  conditions 
even  If  the  system  has  hardening  spring-type  nonlinearity-  One  such 
example  is  shown  In  Figure  12. 


2.4.2  Exponential  Envelope  Function 

For  the  exponential  envelope  function  expressed  by  Eq,  (2-5), 
Eq.  (2-58)  Is  of  the  form 


E[yZ3 


>t  t 

< 4 

sirfcj(j(t-T)s1nwd(t-T)cosB(T-T)  x 

■o  Jd 


n n 

i 

1«1 


n iii 

Z exp[(cw -c4  )t+(cw»-c4 )x-a| t-t}  ]dTdr 
j«1  n 1 n 3 


K e'2cV  P M M 

-2 — z l A,A,  i 

•S  L1  j J1 


I P 

exp[(cwn-a-Cj)T]s1nwd(t-T)  x 


• i 

exp[{cu)n+a“C^)r]s1na){j(t-T)cosB(T-T)dT 


}>  dr  + 


^0 


M M 

z z a4a. 
1 j 1 3 

rt 

■ i 

exp[{cwn+o-Cj)T]s1nu(j(t-T)  x 

0 

( s 

1 

"" 

I 1 

exp(cwn-a-c^)T  sinw(J(t-T)cosB(T-T)dT 

. 

■ (il 

or 

(2-66) 


After  some  tedious  algebra,  we  have 
«,  K_e’2^wnt  K M 


8 


E[yc] 


x Z Z A.A^  E ^ilk^lik 

ul  1-1  j»1  1 3 k«1  13K  13K 


(2-S7) 


where  and  are  defined  In  Eq.  (2-61).  However,  p^,  p2,  and 
p3  are  replaced  by  the  foil  wing. 


P11  * &Va”c1 

plj  " Vcj 

p2i  - cu)n+a^ci 

"31  * "11+"21 
"3J  ■ pij+p2J 


(2-68) 
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Also  all  other  Qj,....Qg,  T.j,....Tg,  Rj,....Rg  are  denoted  by 

Equations  (2-67)  and  (2-38)  nay  be  solved  simultaneously  for 

E[y2]. 

As  a numerical  example,  the  simplest  case  N»1  In  Eq.  (2-6),  that 
1$,  A(t)«e”ct  Is  shown.  In  Figures  9 through  11,  the  normalized  mean- 
square  value  Is  plotted  for  various  value  of  a.  Here,  a,  6 and 
normalization  factor  CQ  are  the  same  as  those  used  in  the  previous 
Chapter  2.4.1.  The  exponential  decay  constant  c of  the  envelope  function 
is  taken  as  c « W. 
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CHAPTER  III 


MUlTID£fiRE£-GF-FREEDOH  NONLINEAR  SYSTEM 


3.1  Statement  of  the  Problem 


Now*  consider  the  H-degnee-of-freedo«  system  governed  by 


y^c^y^fphi  l * ^(t)  1*1 , 2,..., N 


(3-1) 


Or,  this  may  be  written  as 


'ftVi’f ''i1*1  (3‘2! 

where  V is  the  total  potential  energy  per  unit  mass  expressed  by 


* i 2 E E V&frl 

k»l  ’ k*1  ,1*1  * K ,rj 


2 2 2 2 


(3-3) 


The  forcing  function  f ^ ( t)  is  assumed  to  be  represented  by  a product 
of  a well-defined  envelope  function  A^(t)  and  the  uncorrelated  stationary 
random  process  n^(t)  with  zero  mean,  that  Is, 


(3-4) 

(3-5) 

(3-6) 


*j(t)  * Aj{t)n.j{t) 

ECni(t)l  * 0 

E[n1(t)nj(t)]  * 0 if* 

Furthermore  it  is  assumed  that  n^(t)  has  a power  spectral  density  ^(<u) 
which  is  a smooth  function  of  w,  having  no  sharp  peaks.  The  stationary 
response  of  the  system  of  Eq.  (3-1)  to  uncorrleated  white  random  process* 


has  been  studied  by  Caugbey  [7].  He  are  to  determine  the  approximate  meai 
.2/ 


square  response  Ey^(t)]  subject  to  the  assumptions  given  in  Eqs.  (3-4)  th 
(3-6)  by  equivalent  linearization  technique.  The  envelope  functions  A<(t 


rri’W'VWsr't^' 
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considered  In  this  Chapter  are  the  unit  step  function  and  the  exponentla' 


function,  i.e.. 

A1Ct)  « u(t) 

(3-7) 

A (t)  - A^'S* 

(3-8) 

3.2  Response  Formulation 

Let  us  assume  that  an  approximate  solution  to  Eq.  (3-1)  can  be 
obtained  from  the  linearized  equation 


yi+281e^1+w1ey1  * C3"9) 

2 

where  $ Is  the  equivalent  linear  daaplng  coefficient  and  Is  the 
equivalent  linear  stiffness.  Then  the  error  caused  by  this  llnearizatloi 
Is  obtained  from 

•j(t)  « 2(c1wi-B1e)y1+(wJ-wJe)y1-u^-  (3-10) 


where 


n n 


u * k E 


2 2..2„2 


b l a>4ui*y iV^ 

j*1  k*l  j 

Hininizlng  E[e2]  with  respect  to  and  a>1e,  we  obtain 


S^lEjj^^ylHCy^lEEy^ayl] 
ECyfrt^HECy,*,))2 


3U/. 


281e  * 25^^+p 


(3-11) 


(3-12) 


*c 


2 , Eri^E[y,  W/3y<1-Etyty1lE[y13u/3y<] 
1 + E[yf]E[yf]-(E[y1j1])2' 


(3-13) 


It  Is  easy  to  show  that  the  conditions  given  by  Eqs.  (3-12)  and  (3-13) 

2 

vleld  the  true  minimum  of  Ef>£]  If  we  apply  the  sane  argument  used  In 
Chapter  2.2. 

In  order  to  express  the  right  hand  side  of  Eqs.  (3-12)  and  (3-13) 
In  terms  of  the  mean-square  of  the  displacement  and  the  velocity,  we 
must  know  the  joint  probability  density  function  p{yj,  y2,  ..,yN»  y-j, 
y2, — yN).  Since  the  Inputs  f^(t)  are  Gaussian  and  the  nonlinearltles 
of  the  system  are  assumed  to  be  small,  the  outputs  of  the  linearized 
system  are  also  assumed  to  be  Gaussian.  The  displacement  and  the 
velocity  of  the  1-th  mode  are: 


h^(t-T)f^(t)dt 


(3-14) 


[28] 


Mext,  let  us  find  the  covariances  E[yj(t)y.|(t)]  and  E[yj(t)yj(t)]. 


E[y^(t5yj{t)] 


t 

) 


(t~T)hj(t"T)E[f^ (f)f  j(t)3dTdr 


(3-20) 


since  fj(t)  and  f^(t)  are  uncorrelated,  for  i f j,  l.e., 

ECfi(t)fJ(t)3  « A^(t)Aj(t)E[n.j{t)nj(t)]  « 0 i f j (3-21) 
It  Is  concluded  that 


Ety^tJyjft)]  * 0 1 t j 

Similarly,  we  can  show  that 


(3-22) 


EtyWfyt)] 


i 

6 


h^(t-x) 


Shj(t-T)  s i 

" al E[f j(T)fj(t)]dtdr  * 0 

(3-23) 


The  displacements  and  velocities  between  different  nodes  are  Ritually 


^correlated.  Therefore,  the  covariance  matrix  becomes: 


K* 


ECyf]  ELy1y1]  o 

Ety1y13  ELyfl  o 

o 0 ECyjh 

0 0 ECy^] 


0 

0 


» 

£[$  0 


0 

0 

0 

0 


0 0 
0 0 


0 0 
0 0 


E[yjj]  ECy?fyH3 

ECy^]  E[yJ] 


(3-24) 
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i 

\ 


t* 


For  the  covariance  matrix  expressed  by  Eq.  (3**24),  the  2H-fold 
probability  density  function  is  given  by 

K 

p(y] .y^.-.yjj.y! »y2t**.yN)  Pk(yk»yk)  (3-25) 

where  Pk(yk,yk)  is  the  probability  density  function  associated  with 
the  k-th  node  defined  by 


'^(-Vk^^k-Vk)  (3-Z6> 


•tetd^)  - E[yZ]E[yZ]-(Etykyk]) 


ECyf] 


«L  * 


k T&i WJ 


b 


E[ykyk] 
k " mt(\ 9' 


(3-27) 


c.  * 


k 23ecTK£J 


Using  Eqs.  (3-25)  and  (3-26),  E[y^  —4  and  E[y^  can  be  evaluated  as 

yi  y1 


dU 


Efy  JJL*j 

Ly1  3y,J 


) 


3U  N 


•09  4 

2K-fold 


y1  1y7  kl}  Pk{ykyk)dykdyk 


“ 4E[y1]  “kE[y*]+3w^E[y2]} 

kfl 


* w^Etyf]  { E (^[yj^u^Cy*]} 


(3-28) 


[30] 


r 


Similarly, 

y1  k^pk^k^k^dyk^k 

■00 

2*1- fold 

£ 

- { E o^CyJ]  ♦ 2«fECyf]>  (3-29) 

! 

Substituting  Eqs.  (>>28)  and  (3-29)  into  Eqs.  (3-12)  and  (3-13),  mo 
obtain, 

| 2S1e  « 2c1wi  (3-30) 

|| 

II 

J >4e  * <^Ety|l  + 5SwfECyfl)l  (3-31) 


E«,  #4 


These  results  are  the  sane  as  the  equivalent  linear  damping  and 
stiffness  for  the  stationary  processes. 

In  this  Chapter,  the  mean-square  response  Is  formulated  by  use 
of  Eq.  (2-40).  The  mean-square  response  of  the  i-th  mode  is  given  by. 


E[y?]  * 


1 *,(*) 

t 

0 

0 Jo 

V 

exp  {-Cju.j(2t-T-T)  slnw^U-x)  x 


V ■ I I 

s1nwjd(t-T)A.j  (t)A^  (T)cosw{T-T)dTdtd« 


(3-32) 


who:*  * is  defined  by  Eqs.  (3-17).  Now  by  substitution  of  Eqs.  (3-30) 

2 

and  (3-31)  i ‘0  ~q.  (3-17),  is  expressed  as  follows: 


rtf  M»>.  ihm*wik, 


N 

<4d  * + t»(^E  + ?«^Cy?3)3 


(3-33) 


3,2.1 


Substituting  Eq.  (3-7)  into  Eq.  (3-32)  and  integrating,  Eq.  (3-32) 


reduces  to 


E[y*l  « ♦1(M)|H1(n)|2K1(tt,t)dM 


(3-3«) 


where 


(u,t)  « l+e  11  (1+  -p-  sin2ui(jt)+e  1 sin  c^t 


^jWj-Wjj+O)  -CjfUjt  £*W. 

( 1 1 ^ - >-2e  1 1 (cosu)dt+  JJ.  Slnw1dt)  x 
“id  1(1 

2 ’^“i* 

cosut-  ~ slnwx  .tsinwt 

“id  id 


|H1(o>) 


(w^w  r+f&owj^)* 


(3-35) 


(3-36) 


Since  ^(to)  has  been  assumed  to  be  a smooth  function  of  w,  having  no 
sharp  peaks,  and  if  ^ Is  small,  then  the  Integral  of  Eq.  (3-34)  can 
be  evaluated  approximately  by  [1]J 


E r77T;2~T7  Cl’e  0+  ~T~  sin  “idt+  ^7  sin2t!,1dt):I 


(3-37) 
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If  fy(w)  is  given  and  Eqs.  (3-31)  and  (3-33)  are  substituted  Into 
Eq.  (3-37),  we  have  N simultaneous  nonlinear  algebraic  equations  for 

E[.Vj]  (1  * 1,  2, N).  These  can  be  solved  numerically  by  the 

Newton-Raphson  method. 

Kow  consider  the  particular  case  in  which  n^(t)  is  white  noise. 
If  we  denote  4^(w)  * constant  * X,  Eq.  (3-34)  becomes 


E[yf] 


.2  2 


u -2CjWjt  P 

— [1-e  1 1 (V  — p sin  w 

4?i“i“le  “Id 


id 


t+  sin2wid1)] 


u 


id 


(3-33) 

p 

Let  us  now  show  that  the  transient  mean-square  response  E[y^]  does 

2 

not  exceed  the  stationary  mean-square  response  E[yj]$  if  n^(t)  is  white 
noise.  The  stationary  mean-square  response  Is  obtained  by  letting  t -*■  ®. 


- ~r 

^111e 

Substituting  Eq.  (3-31)  into  Eq.  (3-39),  we  find: 

Etifo  (Hu  £ - tA  ■ 0 

Using  Eq.  (2-47),  we  have  the  following  Inequality: 
E[y2]  « 

4c1w1wie 


(3-39) 


(3-40) 


(3-41) 


Substituting  Eq.  (3-31)  Into  Eq.  (3-41)  and  rearranging  the  terms  we 
obtain 


p 


Y 

> 


y- 


%■*, 


,33] 

N 

E[y?]  U+u(  E A[y^]+2wfE[y^)}  - — ^ 1 0 (3-42) 

1 k=l  * * 1 1 4^ 

After  eliminating  irKM^Wj  from  Eqs.  (3-40)  and  (3-42)  and  after  some 
manipulation, 

N N 

(ELvf]-E[yf]sH1+2wf[  J ^(E[y^E[y^]s)]  + u £ <aj;E[yj;]}  i 

k*l  k*l 

N 

ECvf]s  ( £ ^(E[y^s-E[y^)}  (3-43) 

k-1 

Suppose 

£0?]  i E[yf]s  (3-44) 

Equation  (3-44)  implies 
N 

I ^(ECy^-ECy^s)  i 0 (3-45) 

Then  the  right  hand  side  of  the  inequality  (3-43)  is  negative  so  the 
left  hand  side  must  be  negative,  too. 

(E[yf]-EOf]s)  i o 

This  is  a contradiction  of  Eq.  (3-44). 

Hence, 

Etyf ] 1 E[yf3s  (3-46) 

Thus,  it  has  been  proven  that  If  n^(t)  is  white  noise  and  the  envelope 
function  is  the  unit  step  function,  then  the  transient  mean-square 
response  does  not  exceed  the  stationary  mean-square  response. 
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3.2.2  Exponential  Envelope  Function 

For  the  exponential  envelope  function  A..(t)  = e'ci  , Eq.  (3-32) 
now  becomes,  after  double  integration, 


ECyf(t)]  = 


00 

4>i  (w)  | H.a(uj)  1 2W.  (a),  t)da> 


J 

0 


( 3-47) 


where 

2 2 2 

0/-  4-  r • ~U)  * J+W 

w,  (w.t)  = e~2ci  ‘(MJthUt)  [---jr — ] 

“Id 


-2X~(t)cos  t-2X.(t)  •—  sinwt} 
3 4 wid 

X-|(t)  = e ri^(l+  j~sin2a).(jt) 


X2(t)  = e'^i^in^.jt  (3-48) 


X3(t)  - e"2rit(cosa;i.dt+  sina3idt) 


X^(t)  = e“ritsino)^dt 


ri  " ?iwrci 


iHiA^wH2  * 


^id+ri”a32)2+(2riw)2 
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If  Is  either  the  same  order  of  magnitude  as  or  smaller,  then 
the  integration  of  Eq.  (3-47)  may  be  approximated  by  the  following 
expression: 


ECyf(t)] 


"ld+r1>  -2c, t 

y--» — e 1 

4r1 (“1 d+r1 ) 


x [1-e  ^**1^(1+  sinw^ ^t+2  — g~~  sin^w^t)]  (3-49) 

w1d  uid 

Letting  -*■  0,  Eq.  (3-49)  reduces  to  Eq.  (3-38). 
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CHAPTER  IV 
CONCLUSIONS 

In  Chapter  II,  the  time  varying  mean-square  response  of  a non- 
linear single-degree-of-freedom  mechanical  system  to  nonstationary 
random  excitation  characterized  by  the  product  of  an  envelope  function 
and  a stationary  Gaussian  random  process  has  been  considered.  A unit 
step  envelope  function  and  an  exponential  envelope  function  are  consi- 
dered in  conjunction  with  both  correlated  and  white  noise  with  zero 
mean.  The  nonlinear  governing  equation  was  linearized  by  the  method  of 
equivalent  linearization. 

For  the  nonstationary  process,  it  has  been  shown  that  the  equi- 
valent linear  damping  coefficient  and  the  equivalent  linear  stiffness 
for  the  system  with  nonlinearities  involved  only  in  displacements  or 
only  in  velocities  are  the  same  as  those  for  the  stationary  process. 

The  mean-square  response  depends  upon  the  coefficients  of  the 
system  equation,  the  shape  of  the  envelope  function,  and  the  parameters 
of  the  autocorrelation  of  the  process  n(t).  It  was  proved  that  for 
white  noise  modulated  by  a unit  step  function,  the  transient  mean-square 
response  never  exceeds  the  stationary  response.  However,  the  mean-square 
response  to  correlated  noise  modulated  by  a unit  step  function  may 
exceed  its  stationary  value,  especially  when  the  power  spectral  density 
of  the  process  n(t)  has  a sharp  peak,  and  its  maximum  value  becomes 
several  times  the  stationary  value. 
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It  has  also  been  shown  that  the  mean-square  response  of  the 
system  with  cubic  hardening  spring-type  nonlinearity  may  be  greater 
than  the  corresponding  linear  system  response  under  certain  conditions. 

In  Chapter  III,  the  analysis  has  been  extended  to  the  N-degree-of- 
freedom  nonlinear  system  for  the  case  of  mutually  isncorrelated  noise. 


NOMENCLATURE 


envelope  function 

constants  of  the  exponential  envelope  function 

functions  of  the  correlation  function  E[y  ],  E[yy], 
and  E[y2]  defined  by  Eq.  (2-29). 

decay  coefficients  of  the  exponential  envelope  function. 

normalization  factor 

expected  value  of  [ ]. 

time  varying  mean-square  response 

time  varying  mean-square  response  of  the  linear  system. 

stationary  mean-square  response  of  the  linear  system 

stationary  mean-square  response  of  the  nonlinear  system 

difference  between  a nonlinear  system  and  its  equivalent 
linear  system 

impulse  response  function  or  weighting  function  of  the 
equivalent  linear  system 

transfer  function  of  the  equivalent  linear  system 
determinant  of  the  correlation  matrix 
constant 

modulation  function  due  to  unit  step  function 
input  random  process 
probability  density  function 

Wci 

autocorrelation  function  of  input  noise  n(t) 
time 


*****»■ 
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■unit  step  function 

* modulation  function  due  to  exponential  envelope  function 

■ displacement  response 

« decay  coefficient  of  noise  correlation  function 
« frequency  of  noise  correlation  function 

■ equivalent  linear  damping 

* Dirac  delta  function 

* functions  defined  by  Eq.  (3-48).  k » 1,2,..., 4 

■ constants  defined  by  Eq.  (2-52) 

* coefficients  of  the  nonlinear  tens  of  y(t) 

■ system  damping  coefficient 

■ circular  natural  frequency  of  the  corresponding  linear 
system 

» equivalent  linear  stiffness 


■ power  spectral  density  of  Input  noise  n(t) 

* d(  )/dt 
« product 

■v 

■ summation 

* approximately  equal  to 
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Figure  1*  Mean-square  Response  of  the  Nonlinear  Systems  with  Various 

Nonlinesrities  to  White  Noise  Modulated  a Unit  Step  Function 
System  Damping  5*0.025 
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Figure  2»  Mean-square  Response  of  the  Nonlinear  Systems  with  various 

Nonlinearities  to  White  Noise  Modulated  by  a Unit  Step  Function 
System  Damping  0.05 


Figure  3*  Mean-square  Response  of  the  Nonlinear  Systems  with  Various 
Nonlinearities  to  White  Noise  Modulated  by  a Unit  Step 
Function*  System  Damping  ?*0.10 


Normalised  Mean- square  Response 


Figure  4 i Mean-square  Response  of  the  Nonlinear  System 
to  White  Noise  Modulated  byExponential 
Envelope  Function*  System  Damping  5*0.05 
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Figure  5*  Kean- square  Response  of  the  Nonlinear  Systen 
to  White  Noise  Modulated  by  Exponential 
Envelope  Function.  Syetex  Daaping  S* 0.05 
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Figure  10*  Mean-square  Response  of  the  Nonlinear  System  to  Correlated 
Noise  Modulated  by  an  Exponential  Envelope  Function. 

System  Damping  £=0.05 
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Figure  11*  Mean-square  Response  of  the  Nonlinear  System  to  Correlated 
Noise  Modulated  by  an  Exponential  Envelope  Function. 

System  Damping  £=0.05 
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Figure  12 1 Mean-square  Response  of  the  Systems  to  Correlated 
Noise  modulated  by  a Unit  Step  Function. 

System  Damping  £=0.1 
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APPENDIX 


Computation  of  ECyy21**1]  and  ECy2"*2] 


ECyy2"*1]  * 


1 


2ir[det(K)]j 


,00 


y y2mf  1 exp ( -ay 2+2byy-cy2 ) dy  4y 


1 


2ir[det(K)] 


I 


yexp(-cy2)(~  y)2"1*1^)* 


x ; (2mH}-li8v)  ( a )V--2vd- 

1 v-0  vl  4b* 


where 


and 


(2mH;-1;2v/  * (2aw-l)(2w)(2m-1)...e(2»-2v*-l}  v«1,2,‘ 


(2iH-l;-l;0)  - 1 


Elw21"*1!  ■ I I /?2 

J 2«[d.t(K)]*  r 4 *(« 


m-1 


0) 


a * (ac-b*)/a 


where 


[54] 


Noting  that 


I?  m!  / a \V_  /I  x a 

L TSP^TTSr  <5>  <5  + :?> 


Eq.  (2)  reduces  to 


ECyy2"*1]  • <E[y2])"  E[yy] 

2nm! 


E[y^3  » 


2Tr[det(K)P 


y2wf2exp(-ay2+2byy-cy2 ) dy  dy 


2*[det(K)P 


exp[(-c+  L.)y2](^)2^ 


/r  ; 

F * v-0  v*  4b^ 


(ELy2])^1 


2*  m! 
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